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Abstract 

Radiative magnetic dipole decays of heavy-light vector mesons into pseudoscalar 
mesons V — > Pj are considered within the relativistic quark model. The light 
quark is treated completely relativistically, while for the heavy quark the l/rriQ 
expansion is used. It is found that relativistic effects result in a significant reduction 
of decay rates. Comparison with previous predictions and recent experimental data 
is presented. 
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In this letter we consider radiative magnetic dipole (Ml) transitions of the 
ground state vector (V) heavy-light mesons to the pseudoscalar (P) ones, 
V — > P7 (in quark model notations, l 3 Si — > l 1 ^ + 7). For this purpose, 
we use the relativistic quark model based on the quasipotential approach in 
quantum field theory. Recently, this model has been successfully applied for the 
description of different properties of the heavy-light mesons, such as their mass 
spectra [1] and rare radiative decays [2]. Our analysis showed that the light 
quark in the heavy-light mesons should be treated completely relativistically, 
while for the heavy quark it is useful to apply the expansion in powers of the 
inverse heavy quark mass I/tuq, which considerably simplifies calculations. 
The first and sometimes second order corrections in 1/tuq are also important 
for the heavy-light meson decay description. Analogously, for the radiative 
decay calculations considered here, the I/tuq expansion is carried out up to 



Preprint submitted to Elsevier Preprint 



1 February 2008 



the second order, and the light quark is treated completely relativistically (i. e. 
without the unjustified expansion in inverse powers of the light quark mass). 
It follows from the obtained results that the relativistic effects give substantial 
contributions to the calculated decay rates. 

The radiative V — > P7 decay rate is given by [3] 



T = — IA4I 2 , where uj = — — . 

3vr 1 1 2M V ' 



My and M P are the vector and pseudoscalar meson masses. The matrix ele- 
ment of the magnetic moment is defined by 
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A = P-Q, 



(2) 



A=0 



where (P | J M (0)| V) is the matrix element of the electromagnetic current be- 
tween initial vector (V) and final pseudoscalar (P) meson states with momenta 
Q and P respectively. 

We use the relativistic quark model for the calculation of the matrix element of 
the magnetic moment Al (2). In our model a meson is described by the wave 
function of the bound quark-antiquark state, which satisfies the quasipotential 
equation of the Schrodinger type in the center-of-mass frame [1]: 

(b 2 (M) p 2 \ r d 3 q 

*^ P W 7^%q;M)fM(q), (3) 
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where the relativistic reduced mass is 



M 4 - (m| - mlf 



and b 2 (M) denotes the on-mass-shell relative momentum squared 
2 [M 2 -(m 1 + m 2 ) 2 ][M 2 -(m 1 -m 2 f] 

b W = up • (5) 



Here and M are quark masses and a heavy-light meson mass, respectively. 

The kernel V(p, q; M) in Eq. (3) is the quasipotential operator of the quark- 
antiquark interaction. It is constructed with the help of the off-mass-shell 
scattering amplitude, projected onto the positive energy states. An important 
role in this construction is played by the Lorentz-structure of the confining 
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quark- ant iquark interaction in the meson. In constructing the quasipotential of 
the quark-antiquark interaction we have assumed that the effective interaction 
is the sum of the usual one-gluon exchange term and the mixture of vector 
and scalar linear confining potentials. The quasipotential is then defined by 
[4] 

V(p, q; M) = «i(p)« 2 (-p)V(p, q; M) Ul {q)u 2 (-q), (6) 

with 

V(p, q; M) = ^ s ZV(k) 7 f 72 + Konf(k)nT 2;M + (k), 

where a s is the QCD coupling constant, D^ v is the gluon propagator in the 
Coulomb gauge and k = p — q; 7m and u{p) are the Dirac matrices and spinors. 
The effective long-range vector vertex is given by 

r„(k) = 7m + ^<V*"» fe" = (0,k), (7) 

where k is the Pauli interaction constant characterizing the nonperturbative 
anomalous chromomagnetic moment of quarks. Vector and scalar confining 
potentials in the nonrelativistic limit reduce to 

KonfM = (1 - e)(Ar + B), V c s onf (r) = e(Ar + B), (8) 

reproducing 

Konf(r) = Konf(0 + Konf(0 =Ar + B, (9) 

where e is the mixing coefficient. 

The quasipotential for the heavy quarkonia, expanded in p 2 /m 2 , can be found 
in Ref. [4] and for heavy-light mesons in [I]. All the parameters of our model, 
such as quark masses, parameters of the linear confining potential, mixing co- 
efficient e and anomalous chromomagnetic quark moment k, were fixed from 
the analysis of heavy quarkonium spectra [4] and radiative decays [3]. The 
quark masses mj = 4.88 GeV, m c = 1.55 GeV, m s = 0.50 GeV, m u ^ = 0.33 
GeV and the parameters of the linear potential A = 0.18 GeV 2 and B = —0.30 
GeV have the usual quark model values. In Ref. [5] we have considered the 
expansion of the matrix elements of weak heavy quark currents between pseu- 
doscalar and vector meson ground states up to the second order in inverse 
powers of the heavy quark masses. It has been found that the general struc- 
ture of the leading, first, and second order 1/rriQ corrections in our relativistic 
model is in accord with the predictions of HQET. The heavy quark symme- 
try and QCD impose rigid constraints on the parameters of the long-range 
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Fig. 1. Lowest order vertex function corresponding to Eq. (12). Radiation only 
from one quark is shown. 
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Fig. 2. Vertex function T^> corresponding to Eq. (13). Dashed lines represent the 
interaction operator V in Eq. (6). Bold lines denote the negative-energy part of the 
quark propagator. As on Fig. 1, radiation only from one quark is shown. 

potential in our model. The analysis of the first order corrections [5] fixes the 
value of the Pauli interaction constant k — —1. The same value of k was found 
previously from the fine splitting of heavy quarkonia 3 Pj- states [4]. The value 
of the parameter characterizing the mixing of vector and scalar confining po- 
tentials, e = —1, was found from the comparison of the second order (1/toq) 
corrections in our model [5] with the same order contributions in HQET. This 
value is very close to the one determined from considering radiative decays of 
heavy quarkonia [3], especially the Ml-decays (e. g. the calculated decay rate 
of J/ty — > r/cj can be brought in accord with the experiment only with the 
above value of e) . 

In the quasipotential approach, the matrix element of the electromagnetic 
current between the states of a vector V meson and a pseudoscalar P 
meson has the form [6] 



where T M (p, q) is the two-particle vertex function and ^v,p are the meson wave 
functions projected onto the positive energy states of quarks and boosted to 
the moving reference frame. The contributions to Y come from Figs. 1 and 2. 
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d 3 p d 3 q 



^ PP (p)r M (p, q)* yQ (q), 



(10) 
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The contribution is the consequence of the projection onto the positive- 
energy states. Note that the form of the relativistic corrections resulting from 
the vertex function explicitly depends on the Lorentz structure of the 
(^-interaction. Thus the vertex function is given by 

r,(p,q) = r( 1 )(p,q) + r( 2 )( P ,q) + ---, (11) 

where 

r«(p,q) = e 1 u 1 (p 1 )^u 1 (q 1 )(2nf5(p 2 - q 2 ) + (1 <- 2), (12) 



and 



^ ) (p,^) = e 1 u 1 (p 1 )u 2 (p 2 )\v(p 2 -q 2 ) 



A^i 



7i°V(p 2 - q 2 ) \u 1 (q 1 )u 2 (q 2 ) + (!<-> 2). (13) 



+7lM ei(A;i) + ei(pi) 

Here ei )2 are the quark charges, ki = pi — A; k[ = qi + A; A = P — Q; 

e(p) - (m 7 + 7 ( 7 p)) 



AH(p) = 



Hp) 



e(p) 



pz _|_ 



It is important to note that the wave functions entering the current matrix 
element (10) cannot be both in the rest frame. In the initial V meson rest 
frame, the final P meson is moving with the recoil momentum A. The wave 
function of the moving P meson fp A is connected with the wave function in 
the rest frame ^ P0 = ^ P by the transformation [6] 

*pa( P ) = D^iRYjD^iRYj^Poip), (14) 



where R is the Wigner rotation, La is the Lorentz boost from the rest 
frame to a moving one, and D 1 ^ 2 (R) is the rotation matrix in the spinor 
representation. 

We substitute the vertex functions and given by Eqs. (12) and (13) 
in the decay matrix element (10) and take into account the wave function 
transformation (14). To simplify calculations we note that the mass of the 
heavy-light mesons M Vt p is large (due to presence of the heavy quark M VjP ~ 
itlq) and carry out the expansion in inverse powers of this mass up to the 
second order. Then we calculate the matrix element of the magnetic moment 
operator (2) and get 
(a) for the vector potential 
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ei(p)+mi e 2 (p)+m 2 



.* v (p) + (l<-2), 
(15) 



(b) for the scalar potential 



= 

J (2tt) 3 



P(PJ 



ei 



2ei(p) 
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.* V (P) + (1^2). 

(16) 



Note that the last terms in Eqs. (15), (16) result from the wave function 
transformation (14) from the moving reference frame to the rest one. It is 
easy to see that in the limit p/m — > the usual nonrelativistic expression for 
the magnetic moment follows. 

Since we are interested in radiative transitions of the vector mesons to the 
pseudoscalar mesons it is possible to evaluate spin matrix elements using the 
relation (o^) = — (cr 2 ). Then assuming one quark to be light q and the other 
one Q to be heavy and further expanding Eqs. (15), (16) in the inverse powers 
of the heavy quark mass tuq up to the second order corrections to the leading 
contribution we get 

(a) for the purely vector potential (e = 0) 



M v = 



m n 



1 1 

+ 



2m q [\e q (p)/ 3 \e q (p)[e q (p) + m q ] \e q (p) M v 



[1 + K)m q 



, e a(p) Vq(p)+ m q m Q , 



2m 



x 2(p 2 ) | l + /t/p 2 



6M v m Q \e q {p) i 
2 \\ 



Q 



3mf 



+ 



3 \m Q \m Q e q (p) + m q 



QM V \m Q e q (p) + m, 



(17) 
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(b) for the purely scalar potential (e — 1) 




m q (M v - m Q ) \ m g / p 2 
e 2 (p) / 6M v m Q \e q (p) 
2 2e g (p) 



\mQ 3[e g (p)+m g ] 3 My [e,(p) + m 9 



(P 2 ) ( 1 + 1 




6mQ \m Q M V/ 



3 My \e 9 (p) + m 



(1* 



Here (■ ■ •) denotes the matrix element between radial meson wave functions. 
For these matrix element calculations we use the wave functions of heavy- 
light mesons obtained in Ref. [1]. It is important to note that in this reference 
while calculating the heavy-light meson mass spectra only the heavy quark was 
treated using the 1 / rriQ expansion but the light quark was treated completely 
relativistically. 

The values of decay rates of mesons with open flavour calculated on the basis 
of Eqs. (1), (17), (18) are displayed in Table 1. In the second column (r NR ) we 
give predictions for decay rates obtained in the nonrelativistic approximation 
(p/m — > 0) for both heavy and light quarks. In the third (r y ) and fourth (r ) 
columns we show the results obtained for the purely vector and scalar confining 
potentials, respectively. And in the last column (r) we present predictions for 
the mixture of vector and scalar confining potentials (8) with the mixing 
parameter e = —1. As seen from this Table relativistic effects significantly 
influence the predictions. Their inclusion results in a significant reduction of 
decay rates (r NR /T = 2 -=-4.5). Both relativistic corrections to the heavy quark 
and the relativistic treatment of the light quark play an important role. The 
dominant decay modes of D* mesons are the strong decay D* — > Dtt, which 
is considerably suppressed by the phase space, and the electromagnetic decay 
D* — > D7. The corresponding branching ratios are known already for a long 
time and listed in PDG tables [7]. However, the total decay rates of D* mesons 
were not measured until recently. In Ref. [8] CLEO collaboration reported 
the first measurement of the D* + decay width T(D* + ) = 96 ± 4 ± 22 keV. 
Combining this value with the measured BR(D* + — > -D + 7) = (1.6±0.4)% [7], 
the following experimental value of the decay rate can be obtained: T(D* + — > 
D + ^i) = (1.5 ± 0.6) keV. Our model prediction is in agreement with this 
experimental value. However, the experimental errors are still large in order to 
discriminate the relativistic and nonrelativistic results. In the case of B mesons 
the pion emission is kinematically forbidden, so the dominant decay mode is 
electromagnetic decay B* — > B^y. None of B* widths has been measured yet. 
We also present our predictions for the radiative decay rate of B* meson, 
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Table 1 

Radiative decay rates of mesons with an open flavour (in keV). 



Decay 




T v 




r 


D*± _ 


-> Z) ± 7 


2.08 


0.60 


0.28 


1.04 


D*° - 




37.0 


14.3 


17.4 


11.5 


D* - 

^ s 


+ D a <y 


0.36 


0.13 


0.08 


0.19 


B*± - 


-» £±7 


0.89 


0.24 


0.29 


0.19 


B*° - 


^B°7 


0.27 


0.087 


0.101 


0.070 


B*- 


^B sl 


0.132 


0.064 


0.074 


0.054 


B*- 




0.073 


0.048 


0.066 


0.033 



which consists of two heavy quarks (b and c). Therefore the expressions (17) 
and (18) can be further expanded in inverse powers of both quark masses up 
to the second order. 

In Table 2 we compare our predictions for radiative decay rates of vector 
heavy-light mesons with other theoretical results. We show the predictions 
obtained in quark models [9-11], in the framework of heavy quark effective 
theory (HQET) combined with vector meson dominance (VMD) hypothesis 
[12] and in QCD sum rules [13-15]. These predictions vary quite significantly 
from each other. Our predictions are in rough agreement with the quark model 
calculations of Ref. [10], with HQET+VMD results of Ref. [12] and with some 
of the predictions of the QCD sum rules. Q It is important to note that in 
our calculations we do not need to introduce the anomalous electromagnetic 
moment of the light quark as it is done in Ref. [11], where it was found that 
its value should be rather large (~ 0.5) in order to get agreement with the 
experimental (CLEO) value for D* + — > D +r y decay rate. \f\ The large value 
of the anomalous electromagnetic moment is not justified phenomenologically 
(see e. g. Refs. [16,3]). The other differences of our calculations from those of 
Ref. [11] are the Lorentz structure of the confining potential and a more com- 
prehensive account of relativistic effects. In particular, the relativistic trans- 
formation of the meson wave function from the rest frame to a moving one 
given by Eq. (14) is missing in Ref. [11]. 

In Table 3 we present the comparison of our results with the predictions of 
different quark models [17-19] for the rates of the radiative Ml-transitions 
(l 3 Si — > l 1 ^ + 7) in the heavy-heavy B c meson. There we also give the pre- 
dicted values of the photon energy, which is determined by the mass splitting 

1 The QCD sum rule results [13] for the D* decays were obtained using the l/m c 
expansion, which could be inaccurate due to the large value of the l/m c corrections. 

2 In Table 2 we show predictions of Ref. [11] for the value of anomalous electro- 
magnetic quark moment equal to zero. 
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Table 2. Comparison of different theoretical predictions for radiative decays of heavy-light mesons (in keV). 

Quark models HQET+VMD QCD sum rules 

Decay our [9] [10] [11] [12] [13] [14] [15] 

£,.±_ >£) ± 7 104 o.36 1.72 0.050 0.51 ±0.18 0.09lg;$ 1.5 0.23 ±0.1 

D*° -> D°7 11.5 17.9 7.18 7.3 16.0 ± 7.5 3.7 ±1.2 14.4 12.9 ± 2 

D* s ^D sl 0.19 0.118 0.101 0.24 ±0.24 0.13 ±0.05 

s *±^£± 7 o 19 26i 0.272 0.084 0.22 ± 0.09 0.10 ± 0.03 0.63 0.38 ± 0.06 

B*° _> S°7 0.070 0.092 0.064 0.037 0.075 ± 0.027 0.04 ± 0.02 0.16 0.13 ±0.03 

5* -> 5 s7 0.054 0.051 0.035 0.22 ± 0.04 



Table 3 

Comparison of theoretical predictions for the radiative B* — > B c -f decay. 





our 


[17] 


[18] 


[19] 


Photon Energy (MeV) 


61 


72 


64 


55 


T(B* -> B cl ) (eV) 


33 


135 


60 


59 



of the vector and pseudoscalar ground states. In previous calculations [17-19] 
the nonrelativistic expression for the matrix element of the magnetic moment 
was used. We see that even in the heavy-heavy B c meson inclusion of the 
relativistic effects results in considerable reduction of the radiative Ml-decay 
rate. As can be seen from Table 1, this reduction is evoked by significant con- 
tributions of relativistic effects for the c quark, since it is not heavy enough, 
as well as by the special choice of the mixture of vector and scalar confining 
potentials in our model (8). 

In summary we calculated radiative Ml-decay rates of mesons with open 
flavour in the framework of the relativistic quark model. In our analysis the 
light quark was treated relativistically, while for the heavy quark the l/m,Q 
expansion was carried out up to the second order. Relativistic consideration 
of the light quark, relativistic heavy quark corrections as well as Lorentz- 
structure of the confining potential considerably influence the predictions. We 
find that only the mixture of vector and scalar confining potentials (8), with 
the mixing coefficient fixed previously from quarkonium radiative decays [3] 
and weak decays of heavy-light mesons [5] , is in agreement with recent CLEO 
data for the D* + — > D +r y decay rate. More precise measurement of this decay 
rate and the measurement of radiative Ml-decays of other heavy-light mesons 
will be crucial for testing the relativistic quark dynamics. 
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